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Magnetic field in Cepheus A as deduced from OH maser 
polarimetric observations 
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ABSTRACT 

We present the results of MERLIN polarization mapping of OH masers at 1665 and 
1667 MHz towards the Cepheus A star-forming region. The maser emission is spread 
over a region of 6 arcsec by 10 arcsec, twice the extent previously detected. In contrast 
to the 22-GHz water masers, the OH masers associated with Hll regions show neither 
clear velocity gradients nor regular structures. We identified ten Zeeman pairs which 
imply a magnetic field strength along the line-of-sight from —17.3 to +12.7mG. The 
magnetic field is organised on the arcsecond scale, pointing towards us in the west 
and away from us in the east side. The linearly polarized components, detected for 
the first time, show regularities in the polarization position angles depending on their 
position. The electric vectors of OH masers observed towards the outer parts of Hll 
regions are consistent with the interstellar magnetic field orientation, while those seen 
towards the centres of Hll regions are parallel to the radio-jets. A Zeeman quartet 
inside a southern Hll region has now been monitored for 25 years; we confirm that the 
magnetic field decays monotonically over that period. 

Key words: masers — polarization — magnetic fields — stars: formation — ISM: 
individual: Cep A — molecules — radio lines: ISM — Hll regions 



1 INTRODUCTION 

The connection between interstellar magnetic fields and the 
phenomenon of bipolar outflow from young stars is a topic 
of great observational and theoreti c al interest (e.g.|Ba chilled 
119961: IWard- Thompson et alJl2000t iKudoh fc Shibatalll997t 
and references therein). OH masers provide a useful tool for 
probing the kinematics and the magnetic field on subarcsec- 
ond scales, ne ar the centres of bipolar outflows from massive 
young stars JCohenlll989T) . The OH masers are generally 
found at distances from a few hundred to a few thousand 
au from the central source, and are sensitive to magnetic 
fields as small as 1 mG through the Zeeman effect. A se- 
ries of OH maser studies has revealed a common pattern 
of OH masers tracing a dense molecular disc perpendicu- 
lar to the bipolar outflow, with the magnetic field twisted 
in the way predicted by the magnetohydromagnetic models 
by lUchida fc Shibatal il985T) and others (Hutawarakorn fc 
Cohen 1999, 2003, 2005: iHutawarakorn. Cohen fc Brebneil 
120021: iGrav. Hutawarakorn fc Cohenl l2003t) . The magnetic 
field is oriented along the outflow direction on the 10-arcsec 
scales probed by submm-polarimetry, but has a toroidal 
component wound around the outflow on the arcsec scale 
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probed by the OH masers. This appears in the form of a 
magnetic field reversal on opposite sides of the disc. In the 
present paper we consider the OH masers associated with 
the well known outflow source Cepheus A. 

Cepheus A is a nearby star- forming region (within 
1 kpc) and the densest part of the molecular cloud com- 
plex located to th e south of the Cepheus OB3 associa- 
tion dSargentlll977f) . The large-scale (~10 arcmin) morphol- 
ogy of this cloud seen in CO emis sion implies bipolar out- 
flow at velocities of ~25kms _1 iRodriguez. Ho fc Morar] 
1980). The base region of Cep A observed in the NH3 
and CO lines with angular resolutions of 2 arcsec and 
15 arcsec respectively, s hows a quadrupolar s tructure of 
the molecular outflow jTorrelles et alJ ^393). The ob- 
servations of several therm al lines confirm the pres- 
ence of multiple outflows IINaravanan fc Walked Il996l : 
ICodella et alJ 120031) . The complexity of the spatial and 
velocity structure of the Cep A region, characteristic of 
turbulent flow, is supported by maps of the H2 emis - 
sion line at 2.1/xm dHiriart. Salas fc Cruz-Gonzalezl feoon. 
iHughes fc Wouterlootl jl984) found several compact Hll re- 
gions (CHIIRs) at centimetre wavele ngths; these mostly fo l- 
low the edges of the NH 3 clouds (iTorrelles et alJ I1993T) . 
All of those CHIIRs are younger than 10 3 yr and some 
of them are likely associated with high-mass powering 
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stars of spectral types B3 or earlier llHughesl l200ll) . 
Some CHIIRs are accomp anied by clusters of OH , H2O 
and CH3OH masers fe.g. ICohen. Rowland fc Blaiil Il984t 
ITorrelles et~aLlll996t iMinier. Booth fc Conwavl 120001) . Sys- 
tematic studies of OH masers with high angular resolu- 
tion revealed that most OH components move away from 
the central CHIIRs on an estimated outflow time-scale 
of ~300yr (iMigenes. Cohen fc Bre bncr 1992). Observations 
of H2O masers indicate the presence of remarkable lin- 
ear and arcuat e structures, which can delineate shock 
fronts or discs (ITorrelles et al ]|l99d: ITorrelles et alJl200lt 
ICuriel et afl 120021 : iGallimore et al]l2003l) . The H?Q maser 
structure centred on the source HW2 is interpreted as due 
to a molecular disc, 600 au in diameter, which is perpen- 
dicular to the thermal radio jet. The inferred mass of the 
powering source is ^ 20 Mq, i.e. it is a BO. 5 or earlier star 
ilTorrelles et al.lll99fj) . 

In the paper we report full polarimetric observations of 
CepA in the 1.6 GHz OH maser lines. Zeeman splitting of 
the OH maser lines at 1665- and 1667 MHz gave a value of 
+3.5 mG, where the plus sign indicates a field directed away 
from us. Another Zeeman pair has been identified in the 22- 
GHz water maser li ne; the inferred field strength is —3.2 mG 
ISarma et aljEool . that is, pointing towards us. Such mag- 
netic f ields are believed to have an non-negligible dynamical 
effect iMigenes et alJll992l) . We attempt to extend previous 
the observations by using new polarization data to obtain 
an overall structure of magnetic field in the central region 
of Cep A. 



2 OBSERVATIONS AND DATA REDUCTION 

CepA was observed in 1999 in all four ground-state OH 
transitions using 6 telescopes of the Multi-Element Radio 
Linked Interferometer Network (MERLIN) (Table 1). Each 
transition was observed in left- (LHC) and right-circular 
polarization (RHC) simultaneously and correlated to ob- 
tain all four Stokes parameters. The spectral bandwidth of 
250 kHz, covering a velocity range of 45kms~ 1 , was divided 
into 128 channels, yielding a channel spacing of 0.35 kms -1 . 
The band centre was set to a local standard of rest (LSR) 
velocity of —14 kms -1 . 

The phase-referencing technique was applied at all 
epochs with the continuum source 2300+638 used as a 
phase-calibrator. Since the reference source was too weak 
to be observed in narrow-band mode, it was observed 
in wide-band mode (16-MHz) at the appropriate frequen- 
cies (given in Table 1) to achieve the optimum signal-to- 
noise ratio. The cycle-times between CepA and the refer- 
ence were 5.5min+2min at the main line transitions and 
6.5min+2min at the satellite line transitions. In addition 
3C286 and 3C84 were observed for 1-3 hours in all ses- 
sions, at the relevant narrow and wide-band frequencies cor- 
responding to the frequencies given in Table 1 (corrected to 
the appropriate Vlsr for the narrow-band configurations). 
3C286 was used as the primary calibrator for the flux density 
scale and for the polarization angle; 3C84 was used as the 
point source, bandpass calibrator and for the polarization 
leakage corrections. 

D ata reduction was car ried out using standard proce- 
dures JPiamond et alJl2003ft with the local d-programs at 



Jodrell Bank and the Astronomical Image Processing Sys- 
tem (aips). To map the emission we used a circular Gaus- 
sian beam of full width half maximum (FWHM) 120 mas 
and the pixel separation was 40 mas in all epochs of obser- 
vations (the synthesised beams for all data are given in Table 
1). The rms noise (a) levels in emission- free Stokes / single 
channel maps were typically a few mjybeam -1 (Table 1). 
The 1665-MHz data taken on May 20 were affected by the 
ringing effect. This data set was therefore smoothed using a 
Hanning function. 

A region of 20 arcsec by 20 arcsec centred at RA= 
22 h 56 m 18!033, Dec= 62°0l'48'.'347 (J2000) was inspected 
for total intensity emission (I Stokes parameter) stronger 
than 5a. We fitted 2 dimensional Gaussian components to 
each patch of maser emission in each channel to measure its 
position and flux density. The emission was considered to be 
real if present at least in three contiguous channels. We mea- 
sured the flux densities in the images of the other Stokes pa- 
rameters at the position of each I Stokes component, using 
the local AIPS task mfquv. We grouped series of components 
above 5a arising from similar positions in successive chan- 
nels into features and analysed them. More details on data 
reduction procedures and determinatio n of the uncertain- 
ties in the maser positions were given bv lNiezurawska et alJ 
l2004h . 



3 RESULTS 

3.1 Polarization spectra 

Spatially complex OH maser emission was found at 1665 
and 1667MHz in the velocity ranges from —26 to 0.1 kms -1 
and from —16 to — 3kms _1 , respectively. The spectra of 
the LHC and RHC polarizations and the linear polarization 
(P = (Q 2 + U 2 ) - 5 ) are presented in Fig. 1. No significant dif- 
ferences were noted in the spectra at the two epochs spanned 
by 10 days. Our circularly p olarized spect r a diffe r consider- 
ably from those reported in ICohen et al] <ll984l) . However, 
CepA is noted for OH fl are activity and strong variability 
l|Cohen fc Brebnerll 19851) . 

No 1612-MHz emission was detected above a level 
of 15mJy over the whole region searched at both 
epochs. This is consistent with single dish observations by 
ICohen. Brebner fc Potteil <ll990l) who found weak and broad 
1612-MHz emission which was resolved-out by MERLIN. 

The 1720-MHz emission appeared on 1999 May 20 as 
completely circularly polarized features with a peak inten- 
sity of about 1 Jy (F ig. 1) which decreased by a factor of 
two after four weeks llNiezurawska et al]l2004|) . We did not 
detect any linear polarization at that transition within a 
sensitivity limit of 15mJy. 



3.2 Distribution of the maser components 

All OH maser features above 5a level are listed in Table 2 
and plotted in Fig. 2. In total, we detected 49 OH masers at 
1.6 GHz lines, 31 features are new (mark ed with superscript 
n in Table 2) and 18 were seen previously ||C^h ^n_etal .119841 ; 
IMigenes et ail 1 19921; lArgon. Reid fc Mentenl l2000l) . The6 
arcsec by 10 arcsec region contains 25 LHC and 16 RHC 
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Table 1. Details of the MERLIN observations. 







Phase calibrator 


Synthesised beam 






Date 


Line transition 


Frequency 


Flux density 


HPBW 


P.A. 


Time 
on source 


Rms noise 
per channel 




(MHz) 


(MHz) 


(mjy ) 


(mas) 


(°) 


00 


(rnjyb" 1 ) 


20 May 1999 


1612.231 


1612.4 


124 


180x140 


-27 


7.5 


3 




1665.402 


1661.2 


128 


149x120 


-7 


6. 1 


5 




1667.359 


1661.2 


125 


148x119 


-11 


6.4 


1 




1720.530 


1720.8 


130 


154x107 


-24 


7.6 


3 


30 May 1999 


1665.402 


1661.2 


128 


136x123 


2 


5.2 


3 




1667.359 


1661.2 


130 


136x122 


-5 


5.1 


4 


16 June 1999 


1612.231 


1612.4 


121 


165x115 


-33 


3.7 


1 




1720.530 


1720.8 


150 


143x109 


-32 


3.7 


4 



Antennas: Defford, Cambridge, Knockin, Darnhall, Mark II, Tabley 



Table 2. Positions of OH masers in CepA observed on 1999 May 30. The newly detected components are marked with superscript n. 
The numbers in brackets indicate relative errors in positions (x0?0001 for RA, xO'.'OOl for Dec). 



Vlsr 


RA 


Dec 


c 

Sp 


Fcctturc 


Vlsr 


RA 


Dec 


Sp 


Fcciti 


(kms- 1 ) 


(22 h 56 m ) 


^OZ Ul ) 


I Tir K — 1 ^ 




(kms -1 ) 


(22 h 56 m ) 




{•>y a ) 




1665-MHz LHC 








1665-MHz RHC 








-25.2™ 


18? 1642(2) 


4y . DK)oyz ) 


U.zoO 


A 


-18.0™ 


18? 1646(3) 


4y. D / oyz) 


U. 1 /4 




-16.1™ 


1850771(2) 


45'.' 809 (8) 


0.296 


B 


-14.2 


18? 1143(1) 


46'.'316(1) 


5.240 


b 


-16.1™ 


1759472(3) 


49'.'829(3) 


0.779 


C 


-14.2™ 


17?9474(4) 


49'.'839(3) 


0.359 


c 


-16.1 


1851148(1) 


46'.'320(1) 


8.880 


D 


-13.1 


18?0976(2) 


49'.'343(2) 


0.359 


d 


-14.4™ 


1851383(2) 


46'.' 466(1) 


3.335 


E 


-12.3™ 


18? 1464(3) 


47'.'030(3) 


0.301 


c 


-13.9™ 


1851324(1) 


46'.' 494(1) 


16.250 


F 


-11.9™ 


18?2134(15) 


46'.'510(13) 


0.066 


f 


-13. 9 n 


18! 0951(8) 


45'.'964(6) 


0.579 


G 


-11.7™ 


17?9699(6) 


50'.'012(5) 


0.130 


g 


-13.9 


17? 9651(3) 


50'.' 006 (2) 


1.426 


H 


-11.6™ 


18?1352(1) 


46'.' 487(1) 


2.070 


h 


-12.3 


17? 8382(2) 


48'.'675(1) 


0.622 


I 


-11.5 


17?8662(5) 


48'.'812(4) 


0.196 


i 


-11.5 


17?8656(2) 


48'.'813(1) 


2.224 


.1 


-11.2™ 


17?8843(1) 


46'.'508(1) 


0.515 


j 


-11.4™ 


17?6966(7) 


52'.'311(6) 


0.195 


K 


-8.0 


17?9978(5) 


49'.' 678 (4) 


0.131 


k 


-8.7 n 


1757184(6) 


49'.' 999 (5) 


0.216 


L 


-7.0 


17?9084(5) 


50'.'6132(5) 


0.090 


1 


-8.6 


17? 8852(1) 


46'.' 484(1) 


2.630 


M 


-2.8™ 


17?9985(4) 


49'.' 667(3) 


0.135 


in 


-7.9 


1759947(1) 


49'.'682(1) 


5.192 


N 


-0.6™ 


17?9111(11) 


49'.'361(7) 


0.423 


n 


-7.9 n 


1758273(4) 


53'.' 190(3) 


0.471 





-0.6™ 


17?9498(1) 


49'.'843(1) 


15.370 


o 


-7.8 n 


18?3843(12) 


47'.'314(10) 


0.098 


P 


-0.6™ 


17? 7829(3) 


53'.' 363 (3) 


1.097 


1) 


-7.7™ 


1759582(11) 


49? 175(6) 


0.209 


Q 












-6.7 n 


17?9478(2) 


49'.'868(1) 


0.583 


R 












-6.2 


1851163(1) 


49'.'386(1) 


4.590 


S 












-6.2 n 


1759488(4) 


52'.' 906 (3) 


0.374 


T 












-6.1™ 


17?8903(7) 


46'.' 497(5) 


0.156 


U 












-6.1 n 


18?0778(9) 


48'.' 88 1(6) 


0.167 


V 












-4.9 


18?0497(2) 


47'.'552(1) 


0.481 


w 












-2.4™ 


17?8407(2) 


49'.' 898 (2) 


0.293 


Y 












-0.7 


17?9520(7) 


49'.' 848 (5) 


0.109 


z 












1667-MHz LHC 








1667-MHz RHC 








-15.7™ 


18? 3852(4) 


42'.' 737(3) 


0.936 


A 


-14.6 


18?1123(1) 


46'.'306(1) 


3.962 


a 


-15.7 


1851125(1) 


46'.'305(1) 


4.831 


B 


-5.4 


17?8358(2) 


49^912(2) 


0.611 


b 


-3.1™ 


1851096(4) 


46^334(4) 


0.580 


C 












-3.1 


17? 8372(1) 


49^899(1) 


2.907 


D 












1720-MHz LHC 






1720-MHz RHC 










-13.3™ 


1756124(1) 


44'.'533(1) 


1.2354 


A 


-15.4™ 


17?6147(1) 


44'.'517(1) 


1.122 


a 



features in the 1665 MHz line, four LHC and two RHC fea- 
tures in the 1667 MHz line and single LHC and RHC features 
in the 1720 MHz line. 

For the assumed di stance of 725 pc 
(Blaa uw. Hiltner fc Johnson! Il959l) the projected area 
containing OH maser emission is 4350 x 7250 au 2 i.e. about 



twice as great as that reported previously. The masers are 
mostly projected ag ainst HII regions label l ed as 2, 3b, 3c 
(the nomenclature of lHuehes fc Wo uterloot l dl984l)) a nd 3di, 
3dii, 3diii regions resolved bv iTorreUes^tal .1 lll998l). Note 
that w e use the convention presented bv iToreeUeset alJ 
lll998fi in their Table 1 but not in their Fig. 1. 
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Figure 1. The spectra of Cep A taken with MERLIN in May 
1999. OH main-line spectra are averaged from two epochs (20 
and 30 May). The bottom panel shows the 1720-MHz emission 
detected in 20 May 1999. 



24 OH maser features at 1665- and 1667-MHz are cen- 
tred around the HII region 2 within the velocity range from 



-25.2 kms -1 to -0.6kms" 



14 OH features from both 



transitions with velocities from —16.1 kms -1 to —3.1 kms 
were projected towards the continuum sources 3di, 3dii, 
3diii. Three OH 1665-MHz features with velocities —11.2, 
—8.6, —6.1 kms -1 were close the east edge of the HII re- 
gion 3c. Six features, which all were detected for the first 
time, are not obviously associated with any known contin- 
uum source. Four 1665-MHz masers with LSR velocities of 

— 11.4, —7.9, —6.2 and — 0.6kms _1 are located 3.8 arcsec 
northward of the continuum source 2. One 1665-MHz fea- 
ture with an LSR velocity of —7.8 kms -1 lies about 2 arc- 
sec NW of the 3dii centre, while one 1667-MHz feature, at 

— 15.7kms _1 , is placed 3.8 arcsec southward of the source 
3dii (Fig. 2). Two 1720-MHz features are projected towards 
the southern edge of source 3b and coi ncide with the OH 
4765-M Hz emission which we reported in lNiezurawska et alJ 
(2004). 



3.3 Polarization properties and magnetic field 

We found ten Zeeman pairs with features coinciding spa- 
tially within 36 mas. Table 3 summarizes their properties. 
Seven pairs were found at 1665 MHz, two at 1667 MHz and 
one at 1720 MHz. Five of them, named Zi, Z2, Z3, Z4 and 
Zg are located around the HII region 2. Three pairs, Ze, Z7 
(1665 MHz) and Z7 (1667MHz), are projected at the sources 
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Figure 2. The OH maser distribution in the Cep A region mea- 
sured with MERLIN. Symbols trace the 1665-, 1667- and 1720- 
MHz LHC and RHC features as indicated in the top-left cor- 
ner. Each feature is labelled by its central LSR velocity (italic 
numbers). In a case of Zeeman pairs the dema gnetized veloc- 
ity is given. Ellipses trace the continuum sources iTorrelles et alJ 
1998). Upright numbers labe l the continuum sour c es acc ording to 
iHuehes fc Wouterlootl <1984ft and ITorrelles et al.l <199«ft . 



3di and 3diii. Pair Z5 lies close to the region 3c, whereas the 
only one Zeeman pair at 1720 MHz line appears in the vicin- 
ity of the HII region 3b. We calculated the magnetic field 
strength (B) from the velocity separation of the Zeeman 
components ( AV) using coefficients tak en from Table 2 in 
iD^rvjes] (|l974h . However, as we stated in lNiezurawska et alJ 
(2004), in the case of 1720-MHz transition we assumed that 
the emission came from well— separated a +1 and a -1 com- 
ponents and used the line splitting of 118 kms -1 G _1 . 

Two Zeeman pairs, one at 1665 MHz and one 
at 1667 MHz, labelled as Z7, coincide spatially and 
in velocity. In the past th ey w ere reported by 
IWouterloot. Habing fc Herman! dl980l) as a Zeeman 
quartet. Channel maps of that quartet in linearly polarized 
emission are shown in Fig. 3 and the polarization profiles 
in Fig. 4. The emission at —15.8 kms -1 appears at both 
transitions with similar linear polarized flux densities, 
P ~145mJyb _1 , while the polarization angles differ 
significantly by 23°. We think that blending of features 
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Table 3. 1.6-GHz Zeeman pairs found in Cep A. The peak flux densities (S p ), the separation in velocity (AV) between two features of 
opposite circular polarization, the demagnetized velocities (V<j) of those features and calculated magnetic field strengths (B) are given. 
The negative sign corresponds to the field towards the observer. The last column identifies features according to the nomenclature from 
Table f. 



Polzn. 



RA 



Dec 



Vlsr 



AV 



B 



(22 h 56 m ) (62°01') (kms" 1 ) (Jyb" 1 ) (kms" 1 ) (kms" 1 ) (mG) 



1665-MHz 
Zi 



LHC 
RHC 



17?9651 
17?9699 



50" 006 
5tf'012 



-13.9 
-11.7 



1.426 
0.130 



-2.2 



-12.8 



3.7 H 



Z 2 



LHC 
RHC 



17?9478 
17?9498 



49"868 
49^'843 



-6.7 
-0.6 



0.583 
15.370 



-6.1 



-3.7 



10.3 R 
o 



Z 3 



LHC 
RHC 



17?9472 
17?9474 



49^829 
49"839 



-16.1 
-14.2 



0.779 
0.360 



-1.9 



-15.2 



3.2 C 



Z 4 



LHC 
RHC 



18? 1642 
18? 1646 



49"668 
49"673 



-25.5 
-18.0 



0.285 
0.174 



-7.5 



-21.; 



12.7 A 



LHC 
RHC 



17?8852 
17?8843 



46^'484 
46^'508 



-11.2 



2.630 
0.515 



2.6 



-9.9 



-4.4 M 
j 



LHC 
RHC 



18? 1324 
18?1352 



46^'494 
46" 487 



-13.9 
-11.6 



16.250 
2.070 



-2.3 



-12.; 



3.9 F 
h 



LHC 
RHC 



18? 1148 
18? 1143 



46^'320 
46V316 



-16.1 
-14.2 



8.880 
5.240 



-1.9 



-15.2 



3.2 D 
b 



1667-MHz 

LHC 17?8372 

Z 8 RHC 17?8358 

LHC 18?1125 

Z 7 RHC 18? 1123 



49"899 -3.1 2.907 

49'.'912 -5.4 0.611 

46'.' 305 -15.7 4.831 

46^306 -14.6 3.962 



2.3 -4.3 -6.5 D 

b 

-1.1 -15.2 3.1 B 

a 



1720-MHz 

LHC 17?6124 44V533 -13.3 1.050 2.1 -14.4 -17.3 A 

Z 9 RHC 17?6147 44?517 -15.4 1.122 a 



unresolved by MERLIN is the most likely explanation 
for that difference (see Fig. 3 and compare e.g. panels 
Vlsr=— 14.7kms _1 at both transitions). Moreover, the 
single dish spectra of the quartet showed asymmetries in 
the line profiles, i.e. more tha n one Gaussian c omponent 
was needed for an accurate fit llCohen et al.lll990lb 

Table 4 gives the polarization characterictics of all 
maser features detected in 1999; the LSR velocity, the 
J2000 coordinates together with their errors, the I, Q, U, V 
and P flux densities, the polarization position angle \ — 
0.5 x arctan(C//Q) and its error a x , the percentage of lin- 
ear polarization, mi = (Q 2 + U 2 ) ' 5 /I, the percentage of 
circular polarization, m c = V/I and the percentage of to- 
tal polarization, m t = (rn 2 + m 2 ) ' 5 . The random errors a x 
due to thermal noise were calc ulated according to the for - 
mula a x = 0.5 xap/Px 180°/7r iWardle fc Kronber Jl97l . 
The errors arising from uncertainty in the measurement of 
the polarization angle of 3C286 were 4° and 3° at 1665 and 
1667 MHz respectively. These systematic errors do not af- 
fect comparisons within the same data set, so they are not 
included in Table 4. 

We found 14 linearly polarized OH maser features 



at 1665 MHz and two at 1667 MHz above the 5a level 
(P 0.010 mjy). Even when linear polarization is de- 
tected in Cep A the percentage polarization is generally 
low (^16percent); we detected only two 1665-MHz com- 
ponents with mi of 22 and 31 per cent. In contrast, the cir- 
cular polarization percentages are typically higher with a 
median value of 77 per cent and with two cases fully polar- 
ized, in 1665 MHz line. Similary, the 1667-MHz line shows 
^6 per cent linear polarization but almost full circular po- 
larization with a median value m c =90 per cent. The distri- 
bution of linearly polarized components together with their 
polarization characteristic are presented in Fig. 5. We did 
not find any linear polarization in the 1720-MHz transition 
above the 5a limit. 

4 DISCUSSION 

4.1 OH masers and radio continuum sources 

The present study reveals that the OH masers are spread 
over an area of 6 arcsec by 10 arcs ec. In addition to 
two main clusters of maser components llCohen et alJll984h 
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Figure 3. Channel maps of linear polarization of the Zeeman quartet at 1665 and 1667 MHz. Contour levels are —1, 1, 2, 4, ... X 60 mjy b" 1 . 
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Figure 4. Polarization profiles of the OH Zeeman quartet in 
CepA at 1665 and 1667 MHz. In the lowest panel the circle and 
cross symbols correspond to m c and mi, respectively. 

associated with the radio continuum sources 2 and 3 
(Hugh es fc Wouterlootl Il984h we detected several compo- 
nents which have not been reported before and which are not 
accompanied by radio continuum background. 94 per cent of 
the new components are much brighter than the 0.1 Jy sensi- 
tivity limit of previous observations, demonstrating the high 
variability of OH masers. 

OH 1665-MHz masers around region 2 (which has 
a complex structure) encompass a wider range of veloci- 



ties ( — 25. 2 km s 1 to — 0.6kms x ) than that reported by 
ICohen et all dl984h . This maser complex extends 3 arcsec 
(2175 au) in an E-W direction and its geometrical centre 
is at RA=22 h 56 m 17 i :9609, Dec=62°01'49'.'578 (J2000). That 
gives a distance of 147 mas (107 au) from t he geometrical 
centr e of a disc traced by water masers (Torrell es et alJ 
Il996ft . Although the E-W elongations of the OH and the 
water masers are similar, the position— velocity diagram for 
OH masers (after demagnetization of the velocities to com- 
pensate for Zeeman splitting) does not show the regularity 
expected from a disc scenario. Our results imply that the ex- 
isting disc has to be strongly disturbed at the distances 1-2 
arcsec from the centre of region 2 by the expanding mate- 
rial. In fact, there are molecular outfl ows at the edge (abou t 
1 arcsec from the centre) of region 2 iMigenes et alJ ll992). 

The OH 1665- and 1667-MHz masers projected against 
the southern continuum source 3d cover mainly a ve- 
locity range from — 16.1kms _1 to — 11.6kms _ . How- 
ever, we fou nd a highly red-shifte d component of velocity 
— 3.1 kms" 1 . iTorrelles et alJ <|l99cf) resolved 3d into a chain 
of four individual 22 GHz continuum sources. They point to 
3dii as harbouring a young stellar object. We only detected 
OH maser emission from the west side of this source. If this 
was powered by an outflow we would expect to see a velocity 
gradient or a systematic spatial distribution of the masers 
but this was not appa rent even after demag netization of the 
component velocities. ITorrelles et alJ dl99ct) reported a sim- 
ilar situation for the water masers around that region, there 
were neither clear spatial nor velocity trend in distribution 
of maser spots. 



4.2 Magnetic field 

The most important result of our polarimetric observations 
concerns the overall regularity of the magnetic field in 
CepA. Fig. 5 shows the distribution of maser components 
with polarization information. The magnetic field, derived 
from the Zeeman splitting of OH 1665- and 1667-MHz lines, 
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Table 4. The linearly polarized OH 1665-MHz and 1667-MHz maser components in CepA. The errors are given in brackets and are 
xO.0001 sec for RA, xO.001 arcsec for Dec and in degrees for \- The last column identifies features according to the nomenclature from 
Table 1. 



Vlsr RA Dec I Q U V P x m m c mt 
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is negative (pointed towards us) in the west and positive 
(pointe d away from us) i n the east side of the mapped 
region. ISarma et alJ (120021) detected a Zeeman pair in the 
water maser line at 22 GHz and derived an upper limit to 
the field strength of —3.2 mG. That component lies just on 
the west side of CepA and is consistent with the magnetic 
field configuration inferred from OH masers (Fig. 5). Fur- 
thermore, a Zeeman pair found in the 1720-MHz OH line 
iNiezurawska et al.l l2004) in the SW region of the mapped 
area supports this picture of an organised magnetic field on 
arcsecond scales. 

The magnetic field reversal in CepA is clearer than 
the reversals found in the molecular outflow sources stud- 



ied with MERLIN by Hutawarakorn fc Cohen (1999, 2003, 
2005) and iHutawarakorn et alJ l|2002tl . However it is more 
difficult to relate the organised magnetic field structure seen 
in CepA on the 1 arcsec scale to the diverse outflows seen 
on larger scales from ~1 to 3 arcmin. Table 5 lists the main 
molecular structures reported towards Cep A. There are four 
preferred directions: -45°, 0°, +45°and +90°. Only H 2 S 
shows an outflow at 0° in the expected direction orthog- 
on al to the magnetic field r eversal according to the model 
bv lUchida fc Shibatal l|l985l ). The multiplicity of outflow di- 
rections suggests that there may be several sources power- 
ing them. The present OH maser observations suggest that 
CepA contains at least two active centres. There is a clear 
need for sensitive molecular observations at the arcsec scale 
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Figure 5. An overview of the magnetic field in CepA detected using MERLIN. Numbers indicate the magnetic field strength in mG. 
The directions of the planes containing the electric field vectors are marked by bars with lengths proportional to the percentage linear 
polarization of each fea ture. Triangles trac e the OH 1665- and 1667-MHz Zeeman pairs from this paper. Diamonds show the water 
maser Zeeman pair from lSarma et alJ <2002ft and OH 1720-MHz Zee man pair iNiezurawska et all2004) . The contours are at 0.5 mjy b 
(3 per cent of the peak intensity) of the 22-GHz continuum emission iTorrelles et al.lll99Sft . 



which could disentangle the outflows from the complex re- 
gion at the core of Cep A. 

The electric vectors of linearly polarized features show 
a systematic trend. The polarization position angles \ °f 
components observed in the direction of the outer parts of 
the HII regions are significantly rotated relative to those ob- 
served towards the centres (Fig. 5). The flux-averaged \ of 
the six electric vectors from the outer parts is —87° ±3°, 
which implies that the component of the magnetic field 



vector in the plane of the sky has a direction of 3° ±3°. 
That is closely aligned with the interstellar magnetic field 
orientation ( — 10°) derived fro m the infrared polarimetry 
towa rds the regions 2 and 3 llJones. Woodward fc Kellevl 
2004). The flux averaged \ from the centre of northern re- 
gion 2 is — 37°±1°, which gives 53°±l°for the magnetic field 
orientation in that area. T hat is parallel to the radio jet 
seen at 22 GHz (P.A .=44°) llTorrelles et alJll996l) . The flux 
averaged \ of 1665-MHz features in southern regions 3di 
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Table 5. Structures in Cep A regions. 



Line 


Scale 
(arcsec) 


P.A. 

(°) 


Synthesised 
beam (") 


Note 


Reference 


NH 3 


42x30 


-45 


5.4x5.3 




Torrefies et al. (19931 




96x24 


+45 










48x36 


-135 








12 CO 


300 


+90 


9.9x6.3 


HV outflow 


Naravanan & Walker (1996) 




180 


+85 




EHV outflow 


CS 


91 


-45 




core 




HCO+ 


60 


+57 


3.2x2.6 


outflow 


Gomez et al. 1999") 


SiO 


120 


-33 


4.7x3.7 


disc 




H 2 S 


85 


+135 


14 


HV outflow 


Codella et al. (20031 




171 


+180 




HV outflow 






~60 


+45 




LH outflow 




SO2 


114x29 


+45 


24 







and 3diii (labels as in Fig. 2) are 78°±2°and -4f°±2°, re- 
spectively. Those imply average magnetic field directions of 
-12° ±2° and 49+2°, respectively. 



4-2.1 Physical conditions 

The magnetic field and the gas densit y relationship in C ep A 
region was established to be B~n 0,4 ijOarav e t ajj|Uj9 6|). At 
the ed ge of NH3 molecular clouds in Cep A, iGarav et alJ 
( 1996) found a magnetic field strength of 0.3 mG and a num- 
ber density of 2xl0 4 cm~ 3 . In this paper, we find that the 
magnetic field strength deduced from OH Zeeman splitting 
has a magnitude of 3.2 to 17.3 mG. This leads to derived 
gas densities in the range from 8xl0 6 to 5xl0 8 cm . Such 
numbers are typical for OH masers in star- forming regions, 
according to the model of lGrav. Doel fc Fieldl il99ll) . Addi- 
tionally, this model predicts that in the case of accelerative 
fields with velocity shifts of 2-3 kms" 1 the 1665-MHz line 
dominates at all number densities. In fact, the masing region 
in Ce p A expands with a velocity of 2.5 kms" 1 l|Cohen et alJ 
1990) and the 1665-MHz emission does appear at least twice 
as s trong as the 1667-MHz line (Fi g. 1). 

iHutawarakorn fc Cohenl (120031) suggested that the evo- 
lutionary path of star-forming regions is seen in their polar- 
ization characteristics. They detected a systematic increase 
in the degree of polarization from the oldest source to the 
youngest one in NGC 7538. In general, Cep A showed highly 
polarized masers with weak linear polarization (typically 
a few percent) but strong circular polarization (typically 
>50 per cent). The median m t for 1665- and 1667-MHz fea- 
tures is 80 per cent. That confirms the young evolutionary 
stage of the Cep A region (similarly to the object IRS 11 
in NGC 7538) . There are no significant differences between 
the polarization properties of regions 2 and 3 what implies 
a similar age for all HII regions. 



4-2.2 Decay 0} magnetic field 

ICohen et, alJ l)l99ffl monitored the Zeeman quartet at 1665- 
and 1667-MHz over a 10-year time span. They present ev- 
idence for a 0.4 per cent yr" 1 decay of the magnetic field 
strength in the OH maser region due to expansion of molec- 
ular gas surrounding a young star. In Fig. 6 we present 



measurements of the velocity separations of the compo- 
nents of the Zeeman quartet over a longer period of 25 
years. We include prev i ously reported data from 1980—90 
l|Wouterloot et alJll98d: ICohen et alJll99d : iFish et al]l2003l) 
and unpublished measurements taken with the Lovell an- 
tenna in the period 1990—95 (Cohen et al., in prep.). 
Our interferometric data gave separations between peaks 
of l.gkms" 1 at 1665MHz and l.lkms" 1 at 1667MHz 
in 1999. We also detected the Zeeman quartet with the 
Nangay antenna in 2002 October and 2004 September 
(Szymczak et al., in prep.). The velocity separations were 
then 1.952, 1.93 kms" 1 (at 1665MHz) and 1.15, l.llkms" 1 
(at 1667 MHz) in both epochs, respectively. Taking the 
whole 25-year period we used least-squares fits to estimate 
the rates of change in line splitting, obtaining —0.0053 ± 
0.0006 kms" 1 yr" 1 for the 1665-MHz line and -0.0013 ± 
0.0007 km s^yr" 1 for the 1667-MHz line. The rates of 
change of field strength implied are —0.009 + 0.001 mGyr" 1 
from the 1665-MHz line and -0.004 ± 0.002 mGyr" 1 from 
the 1667-MHz line. 

The two estimates for the magnetic field decay differ 
significantly of the 2.3cr level (3 percent). This is also ap- 
parent in the lowest panel of Fig. 6. The ratio of line split- 
ting should be 1.67 from theory whereas the observed ra- 
tios are 1.70—1.82. It is most likely those deviations are 
caused by blending (Sect. 3.3). In particular, 1665-MHz fea- 
tures C and c and 1667-MHz feature A all have similar ve- 
locities to the Zeeman quartet. It is imposible to correct 
data for this effect. We take the weighted mean of the 1665- 
and 1667-MHz results as our best estimate of the mag- 
netic field decay. This gives a weighted mean of the mag- 
netic field decay rate of -0.0080 ± 0.0009 mG yr" 1 which is 
-0.24 ± 0.03 per cent yr" 1 . 



5 CONCLUSIONS 

We detected many new 1.6 GHz OH masers in the Cep A 
region, and confirm the strong variability of the OH maser 
emission. OH maser features around HII region 2 have an 
elongated E-W distribution, similar to that of the water 
masers. However, the OH maser kinematics do not show the 
regular position-velocity pattern seen in the water masers. 
The most important result from our full polarimetry 



10 A. Bartkiewicz et al. 



2.1 



<© 1 . 



1.3 



1.2 - 



1.1 



1.9 - 



% 1.8 



1.7 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


1 1 | 1 1 1 1 | 






. . i . . . i . . 


. i 




i i i i i i i i i i i ■ i ■ ■ i ■ i i ■ i i i ■ i i i i 






□ □ . 

A □ x • 

□ n □ □ . 

, , 1 , , . , 1 1 . . 


' 

ft 







1980 1985 1990 1995 2000 2005 
Epoch 

Figure 6. The systematic decrease of the magnetic field in 
CepA shown by Zeeman quartet investigations. The open star 
represents data from this paper taken with MERLIN and dia- 
monds represent data taken with the Nangay telescope in 2002 
and 2004 (Szymczak et al., in prep.). The 1979 data are from 
Wouterloot et al.l il980l) (triangle), the 1982-1990 data are from 
Cohen et alj il990ri (squares), the 1990—1995 data were taken 
with the Lovell antenna (Co hen, in prep . ) (fille d circles). A cross 
indicates the VLA data from lFish et alj 120031) . Bars indicate er- 
rors in measurements. For the details of fitting see text. 



observations is the morphology of magnetic field in close 
surroundings of the young stars. The Zeeman pairs reveal 
a reversal of the magnetic field direction on the arcsecond 
scale, pointed towards us in the west and away from us 
in the east. The electric vectors of linearly polarized fea- 
tures show an additional trend. Those from the outer parts 
of HII regions implied the direction of the magnetic field 
vector closely aligned with the interstellar magnetic field 
orientation derived from infrared polarimetry. In contrast, 
the vectors measured towards the central parts of HII re- 
gions were significantly rotated so as to be parallel to the 
radio jet (region 2). The Zeeman quartet showed a continu- 
ing systematic decrease of the magnetic field strength at a 
decay rate of — 0.24 ± 0.03 per cent yr" 1 . That is almost half 
the rate estimated prevously (ICohen et al.lll990l) . However, 
there are also some effects of blending, which cannot easily 
be removed. Observations with higher angular resolution are 
needed in order to resolve features blended by MERLIN. 
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